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ABSTRACT 

Regenerated high-energy emissions from gamma-ray bursts (GRBs) are studied in detail. If the intrinsic 
primary spectrum extends to the TeV range, these very high-energy photons are absorbed by the cosmic in- 
frared background (CIB). Created high-energy electron-positron pairs up-scatter mainly cosmic microwave 
background (CMB) photons, and secondary photons are generated in the GeV-TeV range. These secondary 
delayed photons may be observed in the near future, which are useful for a consistency check of the intrinsic 
primary spectrum. In this paper, we focus on effects of the CIB on delayed secondary emissions. In particular, 
we show that not only up-scattered CMB photons but also up-scattered CIB ones are important, especially for 
low redshift bursts. They also give us additional information on the CIB, whose photon density is not definitely 
determined so far. 

Subject headings: gamma rays: bursts — infrared: general — radiation mechanisms: non-thermal 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are highly relativistic astrophys- 
ical phenomena located at a cosmological distance. Observed 
gamma-ray spectra are typically highly nonthermal and 
extended to the MeV range or above. The relativistic shock 
scenario is one of the leadin g scenarios to r eproduce suc h 
spectra (see reviews, e.g., (Meszaros 2006i IZhang 20071) 1. 
The internal shock model is one of widely accepted mod- 
els. In this model, the GRB prompt emission is explained 
by electromagnetic radiation from relativistic electrons 
accelerated in shocks generated by collisions among the 
subshells. Theoretically, several emission mechanisms 
of GeV-TeV emission have been proposed. The syn- 
chrotro n self-inverse Compton mechanism (SSC) is one of 
them (Panathanassiou & Meszaro s 19961: iDai & Lu 20021 : 



Peer & Waxman 2004 



Guetta & Granot 20p; 
Casanova et al. 2007D . While such models belong to 
leptonic scenarios, hadronic scenarios are also possible. 
In GRBs, protons can be accelerated up to the ultra high 
energy region. high-ene rgy protons can not only radi- 
ate synchrotron photons dTotani 1998) but also produce 
electron-positron p airs, pions and muons v i a the photo- 
hadronic process dWaxman & Bahcall 1997t I Asano 20051: 
I Asano & Takahara 20031) . Synchrotron radiation by elec- 
trons, pos itrons and muons can a lso contribute to resultin g 
spectra ( Der mer" & Atovan 2004t I Asano & Inoue 20"07h . 
Neutrino detection would be strong evidence of baryon ac- 
celeration a nd expected by future neu trin o detectors such as 
IceCube (I Waxman & Bahcall 1 9971 iGuetta et al. 200 lb 

iDerm er and Atova n2003l : 



Guetta et al. 2004a: 



Murase & Nagataki 2006a; Asano & Nagataki 200' 



Sufficiently high-energy photons, including photons 
originating from protons, make pairs mainly via 77 — > e + e~ 
interaction in the subshells and cannot escape from the source. 
As a result, the intrinsic high-energy cutoff of GRB spectra 
is usually determined by the optical depth for pair creation, 
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which largely depends on bulk Lorentz fact ors of subshells 
dLithwick & Sari 200U iRazzaque et al. 2004 . Even if such 
high-energy photons escape from the subshells, these photons 
may suffer from interactions with the cosmic infrared back- 
ground (CIB) and be largely absorbed especially for high-z 
GRBs such as z > 1 . Hence, the detection of TeV photons 
will be very difficult, unless the GRB location is nearby. The 
secondary electron-positron pairs generated by attenuation 
are very energetic, so that they up-scatter numerous cosmic 
microwave background (CMB) photons by the inverse- 
Compton (IC) process. Such seco ndary photons will be 
observed as de lay ed GeV emissio ns (Chens & Cheng 1996; 
| Dai et al. 2001; IDai & Lu 20021; IGuetta & Granot 2003 : 

~ lAndo 20041: 
The delayed 



Wan g et al. 2004 IRazzaque et al. 2 004; 



ICasanova et"aT~ 2007). 



iFan et al. 2004 
secondary emission is indirect evidence of the intrinsic 
TeV emission as well as a clue to pro bing the intergalac- 
tic magnetic (IGM) field dPlaga 1995b . Such delayed 
secondary emission has been discussed in terms of not 
only the internal shock mode l but als o the external shock 
model dMeszaros & Rees 1994 iMeszaros et al. 1994 
Dermer et al. 20001: IZhang & Meszaros 2001 



Derishev et al. 20021 : IWang et al. 2"00lt IWang et al. 2"004 



Ando 20041) . They could be distinguishable by the distinct 
spectral evolution behavior. 

Photon attenuation due to the CIB is very useful as an indi- 
rect probe of the CIB, whose photon density is not satisfacto- 
rily determined. The direct observation of the CIB is difficult, 
because of the bright foreground emission associated with 
zodiacal light as well as emission from our Galaxy. COBE 
DIRBE and COBE FIRAS have succeeded in highly signif- 
icant detections of the residual diffuse infrared background, 
providing an upper bound on the CIB in the infrared regime. 
On the other hand, galaxy counts have given a lower bound on 
the CIB at wavelengths, where no COBE data are available. 
Despite the dramatic progress in observations achieved by In- 
frared Astronomical Satellite (IRAS), Infrared Space Obser- 
vatory (ISO) and the Submillimeter Common-User Bolomet- 
ric Array (SCUBA), the mid-infrared (MIR) and far-infrared 
(FIR) observations do not reach a level of optical and near- 
infrared (NIR) bands, which can be explained by direct stellar 
emission. In this sense, we have not determined the spectral 



2 



Murase, Asano, and Nagataki 



energy distribution (SED) of the CIB with sufficient accuracy 
yet. For recent reviews, see Hauser & Dwek (2001) and Kash- 
linsky (2006). Stecker et al. (1992) proposed that one can use 
photon attenuation in blazars to determine the intensity of the 
CIB, if we know intrinsic spectra of blazars. Subsequent stud- 
ies used observations of TeV emissions from blazars (for one 
of the latest examples, Aharonian et al. (2006)). Conversely, 
we could obtain information on the intrinsic primary spec- 
trum such as the intrinsic high-energy cutoff, if we know the 
CIB accurately. Similarly to cases of blazars, we can expect 
to make use of GRBs as a probe of the CIB. However, owing 
to uncertainties in GRB intrinsic spectra, the depletion due to 
the CIB in high-energy spectra is hard to be estimated. 

Observationally, GeV photons have been detected from 
some GRBs with the EGRET detector on the Compton 
Gamma-Ray Observatory. Especially, E GRET detected pro- 
longed GeV emissio n from GRB 940217 dHurlev et al. 1994 
and GRB 930131 dSommer et al. 19941) . Although we do 
not know the highest energy in GRB spectra observation- 
ally, theoretical consideration and simple extrapolation of 
GRB spectra enable us to expect TeV photons from some 
GRBs. Furthermore, the tentative detection of an excess 
of TeV photons from GRB 970417a at the 3 a level has 
been claimed with a chance probability ~ 1.5 X 10~ 3 by 
the water Cherenkov detector Milagrito (lAtkins et al. 2000T) . 
although Milagro has not observed such signals so far 
(Milag ro Collaboration 20071) . Another possible TeV detec- 
tion of GRB 971110 has been reported wit h the GRAND 
array at the 2.7er level dPoirier et al. 20031) . Staking of 
data from the TIBET array for a large number of GRB 
time window has led to a n estimate of a ~ la compos- 
ite detection significance dAmenomori et al. 20011) . Fur- 
ther observations of such ver y high-energy gamma -ray sig- 
nals by MILAGRO, MA GIC (iMirzovan et aL 20051) . VERI- 
TAS dHolder et al. 20061). HES S dHinton 20041) and CANGA- 
ROOIII (lEnomoto et al. 2 002) might enable us to detect the 
signals in the near future. However, the photon detection in 
the TeV range can be expected only for nearby events, since 
high-energy gamma-rays will suffer from attenuation by the 
CIB. On the other hand, future detectors such as Gamma- 
Ray Large Area Space Telescope (GLAST) will detect many 
GRBs around the GeV range. They will also enable us to 
discuss not only prompt primary emissions but also delayed 
secondary emissions produced via IC up-scattering. 

It depends on the strength of IGM field whether we can 
observe such delayed gamma-ray signals or not. For electron- 
positron pairs produced by attenuation of 1 TeV primary pho- 
tons, the sufficiently strong IGM field B > 10" 16 G leads to the 
large magnetic deflection of pairs, Ob > 1, and long time de- 
lay of secondary photons, Afg > 10 3 s, so that the secondary 
gamma-ray flux will be suppressed. On the other hand, the 
weak IGM field enables us to have possibilities to detect de- 
layed secondary gamma-ray signals. In this paper, we study 
these two extreme cases, the weak IGM field case and strong 
IGM field case. In the former case, the magnetic deflection 
time is not so important and other time scales such as the 
angular spreading time scale are more important for relevant 
energies. In the latter case, it becomes difficult to observe 
delayed secondary gamma-rays directly from each burst, but 
they contribute to the diffuse gamma-ray background. 

In this paper, we study delayed secondary gamma-ray spec- 
tra from GRBs most quantitatively by numerical simulations 
as well as an approximate formula. We focus on effects of the 



CIB by including the contribution from up-scattered CIB pho- 
tons (hereafter, USIB photons) to delayed secondary spectra, 
which can extend delayed spectra to higher energies but has 
not been emphasized in previous studies. Such a study would 
be important for GRBs that can emit ~ TeV emissions in or- 
der to know the intrinsic feature of the source. In addition, it 
would be useful to obtain information on the CIB more quan- 
titatively. Not only GLAST but also MAGIC and VERITAS 
might detect up-scattered CIB photons in the n ear future. 

This paper is structured as follows. In 32.11 we explain the 
models of intrinsic GRB spectra and we describe the delayed 
emission mechanism in 32.21 The CIB model we use in this 
paper is explained in §2.3. We show the method to estimate 
the diffuse gamma-ray background due to GRBs in §2.4. In 
§3, we show the results. Finally, our summary and discussion 
are described in §4. 

2. THE MODEL 

2.1. Model of Intrinsic Spectra 

Throughout the paper, we focus on long GRBs with the 
typical duration T ~ 10- 100 s. Widths of individual pulses 
vary in a wide range. Typical pulses have the duration with 
6t ~ 0.1 — 10 s and shortest spikes have millisecond or even 
sub-millisecond widths. The internal shock model, in which 
gamma-rays arise by the internal dissipation of relativistic 
jets, can reproduce such wide range variability. However, the 
simple synchrotron model cannot explain several properties of 
prompt emission (see, e.g., Meszaros 2006). In this paper, we 
do not consider these open problems on the prompt emission 
mechanism. The observed photon spectrum is well approx- 
imated by a broken power-law, dN 1 /dE 1 oc (E 1 /E^) for 

E™ <E 1 <E h 1 and dN y /dE 7 oc (E^/E^)' for E* < £ 7 < 
£™ ax , where E™ is the synchrotron self-absorption cutoff, and 
£™ ax is the intrinsic high-energy cutoff. 

The intrinsic high-energy cutoff £™ ax is typically deter- 
mined by the opacity of two-photon annihilation into an 
electron-positron pair. In the internal shock scenario, it is 
easy to see that TeV photons can escape from the subshells 
if an internal collision radius and/or a bulk Lorentz factor are 
large enough dLithwick & Sari 200lh . 

In sufficiently high-energy ranges, gamma-rays due to elec- 
tron SSC, proton synchrotron and charged-meson/muon syn- 
chrotron can contribute to intrinsic spectra via cascade pro- 
cesses. The resulting spectra are complicated and the study 
on them is beyond scope of this paper. We will investigate not 
the intrinsic emission but the delayed emission in detail. The 
delayed emission depends on the amount of attenuated pho- 
tons and would not be so sensitive to the detail of the shape of 
intrinsic spectra with a given £™ ax . Throughout the paper we 

adopt three models with total isotropic energy E lso = 10 53 ergs 
for calculation. Model A: a broken-power law spectrum with 
a = 1 and [3 = 2.2, E* = 300 keV and £™ x = 1 TeV. Model B: 
the same as model A but £'™ ax =10 TeV. Model C: a numeri- 
cally calculated spectrum obtained by Asano & Inoue (2007). 
They perform Monte Carlo simulations including synchrotron 
radiation, Compton scattering, pair creation, synchrotron self- 
absorption and particles originating from protons such as elec- 
trons, positrons, muons and pions. We adopt one of their 
numerical results as model C. For details, see Asano & In- 
oue (2007). Parameters adopted to obtain the spectrum in 
model C are energy per subshell £ s h = 10 5l) ergs, E^ = 300 
keV, an internal collision radius r = 10 15 cm and a Lorentz 
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FIG. 1 . — The intrinsic primary spectra we use in this paper for each model. 
The spectra are normalized by the fluences from a source at z = 1 ■ Model 
parameters are described in the text. 

factor r = 100. The magnetic energy density Ub is assumed 
to be 0.1£/ 7 , where U 1 is the photon energy density in the 
subshell. In Fig. 1, we show the intrinsic spectra for the three 
models we adopt. The second peak of model C in Fig. 1 is 
due to SSC. Above this peak energy, photon absorption due to 
pair creation is crucial. The intrinsic GRB duration (defined 
in the local rest frame) is set to T' = 50 s. 

2.2. Delayed Emission 

For typical GRBs at the source redshift z = 1, most 
high-energy photons above ~ 70 GeV produce electron- 
positron pairs. The produced high-energy pairs cause de- 
layed high-energy photon emission via IC up-scattering of 
CMB and CIB photons. The duration o f such delayed 
emission is det er mined by several ef fects dDai et al. 20021 ; 
Dai & Lu 2001 IWang et al. 2004 IRazzaque et al. 2004 
Ando 2 004; Fa n et al. 2004) ; the angular spreading, IC cool- 
ing and magnetic deflection effects. The angular spreading 
time is expressed as, 



(1+z) 



A 



77 



27 e 2 C 



(1) 



where j e is the Lorentz factor of secondary electrons or 
positrons in the local rest frame in the Robertson- Walker met- 
ric (hereafter, the local rest frame) at each z, and A 77 is the 
photon mean free path. The IC cooling time scale is written 

as, 

he 



At 



ic 



; d+z) 



2je 



(2) 



where fic is the cooling time scale in the local rest frame. If 
the magnetic deflection angle is sufficiently small, the mag- 
netic deflection time is, 



At B : 



■■0-+z)-iic0 2 B 



(3) 



where 9 b = ctic/ri is the magnetic deflection angle and is 
the Larmor radius of electrons or positrons. Note that we have 
implicitly assumed 6b *C 9 ; where 9 j is an opening angle 
of GRB jet. Taking into account the GRB duration T, the 
(secondary) duration time scale is estimated by the maximum 
time scale, Af = max[Af ang , Afic, Ats,T]. Examples of At 
adopting the CIB model of Kneiske et al. (2004, see section 
2.3) are shown in Fig. 2. In cases with a weak magnetic field 
such as B < 10~ (I8-I9) G, the angular spreading time scale is 
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FIG. 2. — The maximum time scales for z = 0.1 and z. = 1 involved in 
calculating the delayed secondary spectra. E e is the electron energy in the 
local rest frame. The IGM field is assumed to be B = 10~ 20 G. 



the most important. Of course, we should note that we treat 
the averaged flux over the duration time scale. 

The IGM field make the situation complicated and more 
careful treatments would be needed to evaluate the accurate 
time-dependent flux of delayed emission . The deleye d emis- 
sion can also be a probe of the IGM field (Pla ga 19951) . but we 
do not focus on this topic in this paper. We will treat the two 
extreme cases. One is the weak IGM field case, where we 
can expect detectable delayed secondary gamma-ray signals 
from GRBs. For demonstration, we use an IGM field with 
B = 10~ 20 G, in which Af# is not so important. The other is 
the strong IGM field case, where it becomes difficult to detect 
delayed signals directly from each GRB. We will consider this 
case later. 

Let us use the local approximation that TeV gamma- 
rays from a GRB interact with the CIB field at the 
same redshift, which is justified by numerical calcula- 
tions later. We can obtain delayed spectra from a burst 
with the source redshift z by an analytic approximate 
formula, whic h is given by ( Blumenthal & Goul d 197C ; 
Dai et al. 20021: IWang et al. 200lE IRazzaque et al. 20041: 
Ando 2004: Fa n et al. 2004 . 



^L_(E 7 ,t,z)= I ' dtJdef d le 



dN^ 



dE^dedtd 



"hc- 



dF e 
dtpdj, 

e -fo/Ar) 

~aT 



-, (4) 



where the electron injection spectrum 



dF e 
dtpdj, 



(E e ,t p ,z) = 2 



dj e 



x(l- 



(5) 



and the photon emission spectrum per unit time due to IC scat- 
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dN~, 2-Krhn e c^ 1 dn r 

-nrhr = ^^ — (e,z) 2yln(2y) 

dxdedtj j e e de 



+ (l+2y)(l-y) + 
_ g 7 (l+z) 

-4e 7< (l-x)' 
_ 4e 7f , 



2(1 +wy) 



(i-y)], 



jn e C' 



2 ' 



(6) 



Here, dF lt i/ dE 7 jdt p is the intrinsic primary gamma-ray spec- 
trum of GRB prompt emission, £ 7 ,■ = 2 7c m f ,c 2 /(l +z) is en- 
ergy of primary photons (where the source redshift is taken 
into account), dnjde is the photon density spectrum of the 
CMB and CIB in the local rest frame, and ro is the classi- 
cal electron radius, t is the given observation time of the 
delayed emission, t p is the time when primary photons are 
released, T is the GRB duration, td is defined by td = t — t p 
and r 77 g (E 7 ,z) is the optical depth against gamma-rays prop- 
agating in the universe. The upper bound of the integra- 
tion over 7 e is determined by the high-energy cutoff of the 



intrinsic primary emission, i.e., j e 



(l+z)E™ x /2. On 



the other hand, the lower bound of the integration over 7c 

is 7 e rain = max[tn e c 2 /2e,((l+z)E 7 /s) 1/2 /2]. We exploit Eq. 

iteratively by substituting (dF 1 / rffrf£ 7 )(exp(T 77 s )- 1) into 
dF 1 j/dt p dE 1 j instead of using the intrinsic primary flux. We 
perform such an iterative method in order to include IC scat- 
tering by generated pairs due to re-absorbed secondary pho- 
tons. 

In the above formula, it is assumed that secondary pairs are 
produced only at the source redshift and cooling time scales 
are evaluated with quantities at the source redshift. The pho- 
ton emission spectrum is evaluated with the initial energy of 
secondary pairs and assumed to be constant while the pairs 
cool. However, pairs are produced at various redshifts and the 
cooling rate becomes lower as they cool. Therefore, we also 
execute numerical simulations including IC scattering and 
pair creation. Based on the pair creation rate at each redshift 
due to the CIB (partially CMB), we follow the time evolution 
of the distribution functions of primary photons / 7 (E 7i ,-), sec- 
ondary pairs / e ( 7e ), and secondary photons fi(E-j) from the 
burst time to the present time. The minimal time (redshift) 
step for / 7 (E 7i ,-) in our simulation is dz = 0.005. The pair 
cooling process is followed with a time step fic/100 until pairs 
become non-relativistic. While / 7 (E 7(J ) decreases monotoni- 
cally with time (or remains constant for lower energy range) 
by attenuation, /a(£ 7 ) does not necessarily change mono- 
tonically especially for high-redshift sources, because of re- 
absorption. IC photon spectra are calculated using the Klein- 
Nishina cross section with the Monte Carlo method used in 
(I Asano 20051) . Our method can precisely treat re-absorption 
of secondary photons and energy loss processes of electron- 
positron pairs. 

2.3. Cosmic Infrared Background 

Gamma-ray absorption due to pair creation in cosmologi- 
cal scales depends on the line-of-sight integral of the evolv- 
ing density of low energy photons in the universe. To 



demonstrate the effect of the CIB on delayed spectra from 
GRBs, we need to exploit some model of the CIB. The 
CIB should be explained by a theory from the first prin- 
ciples, but we are far from this ultimate goal owing to 
poor knowledge about star formation, supernova feedback, 
galaxy merging and so on. So far many models of SED of 
the CIB produced by stellar emiss ion and dust re-radiation 
in galaxies have been constructed dTotani & Takeuchi 20021: 
iKneiske et al. 20021: Istecker et al. 20061) . At lower redshifts 
these models agree with each other basically. At higher red- 
shifts, Stecker et al. (2006) found the larger optical depths 
than previously thought because of intergalactic gamma-ray 
absorption motivated by the recent discovery of active star for- 
mation taking place in young galaxies at high redshifts. Such 
model uncertainties will produce corresponding differences. 

In this paper, we use the CIB model developed by Kneiske 
et al. (2002,2004). They developed the evolving model of the 
infrared-to-ultraviolet metagalactic radiation field, based di- 
rectly on observed emissivities. They specially addressed the 
redshift evolution of the SEDs, which are constructed from 
realistic stellar evolution tracks combined with detailed atmo- 
spheric models (Bru zual & Chariot 19931) . and also taken into 
account effects of re-radiation from dusts and Polycyclic Aro- 
matic Hydrogen molecules in the infrared. Their model pa- 
rameterizes the main observational uncertainties, the redshift 
dependence of the cosmic star formation rate and the fraction 
of UV radiation released from the star forming regions. Here, 
we adopt the "best-fit model" of Kneiske et al. (2004), which 
is consistent with the data obtained from recent galaxy sur- 
veys. In Fig. 3, we show the SED of the CMB+CIB we use 
in this paper. For details, see Kneiske et al. (2002,2004). We 
also assume the ACDM cosmology with VL m = 0.3, J7a = 0.7 
and// = 70kms- 1 Mpc" 1 . 

Given the SEDs, we can calculate the mean free path of 
high energy gamma-rays for pair creation or pair creation rate 
at each redshift. Especially, using the cross section of pair 
creation er 77 , the optical depth of the universe is written as, 



r bkg_ 



= / dz 



cdt 



dz 



1 -cos£ 



dn da 71 
de dcosO 



(/cos t)- — = / </ -1111-1! : (E^,6,e) 

(7) 

For details, see Kneiske et al. (2004). For reference, we plot 
optical depths at z = 0.1 for the other models in Fig. 4. Note 
that the simple power-law fitting formula of Casanova et al. 
(2007) overestimates optical depths in comparison with the 




-2 -1 
log(£ [eV]) 

FIG. 3. — The CMB+CIB radiation field in the comoving frame for various 
redshifts, which we have used in this paper. The data about the CIB field are 
taken from the best-fit model of Kneiske et al. (2004) 
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other models above ^TeV. 

2.4. Diffuse Gamma-Ray Background 

As noted before, if the IGM field is strong enough (B > 
10~ 16 G), it is hard to detect delayed emissions as a source 
connecting with GRB prompt emissions because of large Ats- 
Therefore, such emissions may be detected as the diffuse 
gamma-ray background emission rather than delayed emis- 
sions. Now, we consider this extreme case. 

The observed diffuse gamma-ray emission was found to 
be a power-law in ener gy and is highly isotropic on the sky 
(ISreeku mar et a l. 19981) . but it may not be consistent with a 
simple power-law and The origin of the diffuse gamma-ray 
background from extragalactic sources is also an open ques- 
tion. Blazar is one of the most discussed and promising can- 
didates. Other sources such as fossil radiation from acceler- 
ated cosmic rays during the structure formation might give a 
significant contribution. GRB is one of the brightest astro- 
physical phenomena and also can contribute to the gamma- 
ray background. GRBs as sources for the ^MeV gamma-ray 
background were given by Hartmann et al. (2002). Casanova 
et al. (2007) considered GRBs as sources of the ^GeV-TeV 
gamma-ray background. However, since GRBs are rare phe- 
nomena despite of their brightness, the contribution to the 
gamma-ray background will be small, as is shown later. 

We estimate the diffuse gamma-ray background, indepen- 
dently of the IGM, as follows. The number of GRBs is written 
as, 



Ngrb = 

where the volume factor 
dV c 



dz 



Pgrb(z) dV 
l+z dz' 



4ndf 



dz flo(l +z )Vfl B (l+z) 3 + ^A 
(di is the luminosity distance) and the GRB rate 



Pgrs(z) = Po 



22 + e 3 - 4 * 



(l+z) 



3/2 



(8) 



(9) 



(10) 



Here, we have used the SF2 model of Porciani & Madau 
(2001) for the GRB rate with po = 1 Gpc~ 3 yr _1 , assuming that 
the GRB rate traces the star formation rate in a global sense. 
Guetta et al. (2004) obtained such a value of the GRB rate and 
Liang et al. (2007) also reported a similar value. Even though 
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FIG. 4. — The optical depth of high-energy gamma-rays from a source at 
O.I. 



the actual GRB rate may not be a good tracer of the star forma- 
tion rate dGuetta & Piran 2007l:rLe & Dermer 2007I) . our con- 
clusion about the diffuse background would not be changed so 
much because the main contribution to the background comes 
from bursts that occur at z ~ (1 — 2), the number of which 
is observationally determined. The diffuse gamma-ray back- 
ground due to GRBs is estimated by, 



dF 7 



dz 



dN^ \ dN CRB 



dz 



(11) 



where dN-y jdE-ydA is the observed gamma-ray fluence from 
each burst, which is defined by, 



dN-y 
dE.ydA 



1 dN-y 
And^ dEj 



(12) 



where dN 1 jdE 1 is the photon number spectrum (where the 
source redshift is taken into account) and d p is the proper dis- 
tance to a source. We set z m ; n = and z max = 5. 

We have to note that TeV emissions from GRBs in the in- 
ternal shock model can be expected only in the limited cases. 
In the context of the internal shock model, a sufficiently large 
Lorentz factor and/or large collision radius are required. Al- 
though a fraction of such GRBs that can emit TeV gamma- 
rays is unknown, maybe only a fraction of GRBs are TeV 
emitters. Hence, the contribution of GRBs to the diffuse 
gamma-ray background would give an upper limit, since it 
is evaluated under the assumption that all GRBs have spectra 
extended to TeV energies,. 

3. RESULTS 
3.1. Delayed Gamma-Ray Spectra 

In Figs. 5 and 6, we show total fluences of prompt and 
delayed gamma-rays numerically calculated for various red- 
shifts. In model A (E™' dx = 1 TeV), the maximum energy of 
secondary pairs is at most ~ 500 GeV in the local rest frame. 
Hence, the typical energy of up-scattered CMB photons (here- 
after, USMB photons) is ~ 1 GeV, above which USIB photons 
can make significant contributions to resulting delayed spec- 
tra. Such USIB photons form a relatively flat "slope" in spec- 
tral shape in the ~ (10 - 100) GeV range for z < 3. For z > 3, 
such high-energy slope signature becomes difficult to be seen, 
because secondary photons are absorbed again. 

The above picture is changed, if £™ x is beyond 1 TeV. For 
model B {E™ x =10 TeV), the maximum energy of secondary 
pairs is ~ 5 TeV in the local rest frame. Therefore, the en- 
ergy of USMB photons can reach ~ 100 GeV. Similarly to 
the case of £'™ ax = 1 TeV, the contributions from USIB pho- 
tons are important above ^100 GeV. However, such high- 
energy photons may not reach the Earth, because of dupli- 
cated absorption. The optical depth against high-energy pho- 
tons above ~ 200 GeV exceeds unity for z ~ 0.4. Hence, the 
effect of USIB photons will be buried unless GRBs occur at 
sufficiently low redshifts. 

Although we have discussed CIB effects based on fluence 



in Figs. 5 and 6, flux (F 7 ) of delayed photons rather 
than fluence (£ 7 7 = / dtF 7 ) is often used in discussing de- 
tectability of such photons. Hence, in Figs. 7 and 8, we show 
fluxes that are numerically obtained for a source at z = 0. 1 . For 
reference, approximate results obtained by using Eq. are 
also plotted. The delayed time scale is evaluated using Eqs. 
(Q])-©. In both Figs. 7 and 8, the small difference between 
the numerical results and the approximate ones using Eq. (0]i 
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• The same as Fig. 5, but for model B (£™ x = 10 TeV). 



is seen around the peak of the bump formed by USMB pho- 
tons. This difference will be caused by cooled electrons, as 
explained below. The approximation using Eq. means 
that both pair creation and IC scattering are treated as entirely 
local processes. But high-energy electrons will be produced 
after passing ~ A 77 and IC scattering will occur after passing 
~ Aic = cfic- Thus, these propagating electrons suffer from IC 
losses as well as adiabatic losses and effective electron distri- 
bution would deviate from the expression in Eq. ©• Photon 
emissivity will be also changed correspondingly, and the dif- 
ference between two methods appears. For example, the peak 
of the USMB bump becomes more ambiguous due to such 
losses. Nevertheless, the approach to use Eq. (0]l will usually 
work as a reasonable approximation. 

As is shown in Fig. 7 (£ , ™ ax = 1 TeV), delayed emissions in 
~ (10 - 100) GeV due to USIB photons are prominent. In this 
case, the treatment neglecting the USIB effect is not good and 
leads to underestimation of the delayed gamma-ray flux. Even 
for £™ x = 10 TeV (see Fig. 8), the USIB effect is still remark- 
able. In fact, the contribution from USIB photons is dominant 
above ~ 500 GeV. However, the USIB effect is smaller than 
the case of £™ ax = 1 TeV, since a fraction of delayed secondary 
photons is absorbed again by CIB photons. 

USMB photons form the "bump" shape around the £ 7 ~ 
a few X7 e 2 £cMB/(l +z), resembling the Planck distribution. 
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FIG. 7.— Flux for model A (E* 



1 TeV) obtained numerically and 



obtained using Eq. t[4}. The redshift of a source is set to z = 0.1. The primary 
flux is evaluated for the intrinsic GRB duration T' = 50 s (i.e. time-averaged 
flux over the duration is assumed). On the other hand, the delayed secondary 
flux is evaluated at the observation time t = 10 4 s (i.e. the flux which is 
essentially time-averaged over t = 10 4 s is assumed). 
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FIG. 8. — The same as Fig. 8, but for model B and the observation time 
t = 10 2 s. 

Unless duplicated absorption of secondary photons is signifi- 
cant, the ratio of the fluence due to USIB photons to that due 
to USMB photons should reflect the ratio of the CIB energy 
density to the CMB energy density. This statement is basi- 
cally demonstrated in Fig. 9, from which we can see that 
E*<f> 7 at 100 GeV is about 0.01 times lower than that at 1 
GeV for the normal CIB strength, when the burst occurs at 
z = 0.1. However, for the sufficiently high z bursts, the dupli- 
cated absorption of secondary photons becomes significant, 
where the fluence ratio will largely deviate from the ratio of 
the CIB intensity to CMB one. This implies that we need de- 
layed gamma-ray spectra over wide energy ranges in order to 
see the USIB effect on the delayed spectra, since the cutoff 
of delayed secondary emission depends on the CIB strength 
and the distance to the burst. Note that the number of CMB 
photons increases at high redshifts, while that of CIB photons 
does not change monotonically. Therefore, USMB photons 
are more prominent for bursts at high redshifts. 

In model C, the amount of prompt TeV photons is much 
less than that in model A, and the spectrum of primary prompt 
emission is not expressed as a simple power-law. Neverthe- 
less, the ratio of Ez^ in model C is similar to that in model 
A. Thus, the 100 GeV-1 GeV ratio for low redshift bursts is 
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FIG. 9.— The ratio of El 



at 100 GeV to that at 1 GeV. The thick dot- 



dashed line shows the energy fluence ratios for model A and the thin solid line 
shows those for model C in the normal CIB case, where the best-fit model 
given by Kneiske et al. (2004) is used. The dashed and dotted lines show the 
ratios for model A in the stronger CIB case (5 times as strong as the normal 
case) and weaker CIB case (0.2 times), respectively. 

a good indicator of the CIB strength irrespective of the spec- 
tral shape of prompt emissions, unless the amount of photons 
above TeV is considerable. 

In Figs. 10 and 11, we show time-integrated fluxes, i.e., 
fluences at a given time for the model A from sources at 
z = 0.1 and z = 1, respectively. In both figures, we change 
the overall CIB intensity by a factor 5 and 0.2 in order to 
demonstrate effects of the CIB intensity. As we change the 
strength of the CIB field, the delayed secondary fluence also 
changes correspondingly, which is also shown in Fig. 9. For 
z = 0.1, the cosmic space is essentially optically thin for pair 
creation against CIB photons, so that secondary gamma-rays 
are not completely absorbed. The USIB effect is outstanding 
above ~ (10 - 100) GeV range. Even around the bump, where 
USMB photons are dominant, the fluence changes according 
to the amount of absorbed primary photons, which reflects the 
CIB intensity. We can also see the influence of duplicated ab- 
sorption in the high-energy range, although it is small for the 
case with £™ x = 1 TeV and z = 0.1. 

For z = 1, delayed secondary spectra become more compli- 
cated. Above <~ 70 GeV, gamma-rays cannot reach the Earth 
without attenuation. Hence, delayed secondary gamma-rays 
above this energy are absorbed again and regenerated. If we 
change the strength of the CIB field (as represented in the 
dotted or dashed line in Figs. 10 and 11), the cutoff energy, 
at which the optical depth becomes r^ s = 1, also changes. 
Although the height of the USMB bump is also affected by 
the CIB intensity, this influence will be saturated when pri- 
mary gamma-rays are completely attenuated, like in the cases 
shown in Fig. 11. 

Figs. 12 and 13 show the contrastive results among the 
three models. For low redshift bursts as shown in Fig. 12, 
the high-energy slope feature is produced by the CIB for all 
the models. For high redshift bursts, such high-energy USIB 
slope appears for model A and model C, which is shown in 
Fig. 13. For model C, although a bump signature formed 
by USMB photons is difficult to be seen, the USIB effect is 
still important and delayed secondary photons above <~ 5 GeV 
comes mainly from such USIB photons. On the contrary, in 
the case of model B, USIB photons are almost completely 
absorbed again. Hence, the contribution of USIB photons is 
negligible and it is sufficient to consider USMB photons only 
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FIG. 10. — The fluences of primary and secondary gamma-rays. Model 
A with the redshift z = 0. 1 is used. Fluences of prompt primary (upper three 
thin lines) and delayed secondary emission (lower three thick lines) are time- 
integrated over T' = 50 s and t = 10 4 s, respectively. The CIB strength is 
assumed to be normal (the best-fit model, dot-dashed), stronger (5 times, 
dashed), and weaker (0.2 times, dotted). 
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FIG. 12. — The same as Fig. 10, but for models A, B and C with normal 
CIB strength. The dot-dashed and dotted lines are degenerate for prompt 
primary spectra. 

in this case. 

In Figs. 10-13, we have also shown the short-time sensi- 
tivity curves of GLAST and MAGIC just for comparison. 
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FIG. 13. — The same as Fig. 12, but for z= 1. The dot-dashed and dotted 
lines are degenerate for prompt primary spectra. 

For GLAST, the fluence threshold is roughly proportional to 

1 jl 

(t/2.4 x 10 4 s) for the long-time sensitivity regime (expo- 
sure time t > 2.4 x 10 4 s), and is roughly constant for the 
short-time sensitivity regime (exposure time t < 2.4 x 10 4 s). 
The short-time sensitivity is calculated from the effective area 
of LAT onboard GLAST (which is shown in "http://www- 
glast.slac.stanford.edu/software/IS/glast_lat jperformance.htvri 
under the criterion that at least 5 photons are collected 
(Kamae 2006, Private Communication). Note that the 
GLAST sensitiv it y curves shown in some references 
dPetrv et al. 19991 : iRazzaque et al. 20"04T) are overestimated. 
For MAGIC, the short-time sensitivity is roughly estimated 
from the effective area for the zenith angle with 20 degree 
by the criterion that at least 10 photons are collected, 
although the actual detectabi lity requires careful analyses 
(IMAGIC Collaboration 2006bt) . In addition, the duration 
time of delayed signals can be much longer than that of 
prompt signals. For such long lasting signals, we have 
to take into account the fact that the fluence sensitivity is 
proportional to t '/ 2 for the long-time sensitivity regime. 

For 10 53 ergs bursts at z = 0.1, delayed signals can be de- 
tected by GLAST for model A and model B. In addition, 
the USIB effect around (10- 100) GeV in delayed emissions 
could be detected by MAGIC if the CIB strength is strong 
enough. 

Even for the case of z = 1, MAGIC and GLAST have pos- 
sibilities to detect prompt primary signals (although MAGIC 
has not observed such high-energy emission up to now), if 
GRBs with 10 53 ergs are ~ 10 GeV - 10 TeV emitters. How- 
ever, even GLAST could not see delayed secondary compo- 
nents for z > 1 unless bursts are much more energetic. This re- 
sult is different from that of Razzaque et al. (2004). This dis- 
crepancy comes from just their overestimation of the GLAST 
sensitivity. Only energetic and/or low redshift bursts allow 
GLAST to detect secondary delayed signals from GRBs. 

3.2. Diffuse Gamma-Ray Background 

In the previous section, we have considered the weak IGM 
field case. Here, let us consider the opposite case, where 
secondary gamma-rays contribute to the diffuse gamma-ray 
background. Fig. 14 shows the resulting diffuse background 
in the sense of cumulative gamma-ray background. Even if 
the primary photons are assumed to contribute to the dif- 
fuse background, the contribution is much smaller than the 
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FIG. 14. — The possible diffuse gamma-ray background from GRBs for 
each model. It is assumed that the GRB rate traces the SF2 model with the 
local rate pgrb = 1 Gpc~ 3 yr~' . Note that the contribution from delayed sec- 
ondary emission is expected to be smaller than that from prompt primary one. 
The dot-dashed and dotted lines are degenerate for prompt primary spectra. 

EGRET limit. The contribution due to delayed emissions 
is much less important. This is easily understood as fol- 
lows. If we assume the local GRB rate pgrb (without beam- 
ing correction) ~ 1 Gpc~ 3 yr _1 and the released isotropic en- 
ergy £i so ~ 10 53 ergs, the Hubble time f H (~ 10 10 yr) and 
), the possible cosmological evolution factor on the rate leads to 
the diffuse background ~ (l/47r)c£j so p G RB(z= Ljf H ~ 

10~ 7 GeVcrrf 2 s _1 sr" 1 , which is much smaller than the EGRET 
limit ~ 10~ 6 GeVcm" 2 s" 1 sr" 1 at GeV. Hence, the contribution 
from GRBs to gamma-ray background is expected to be at 
most ~ 10 %. 

4. SUMMARY AND DISCUSSION 

In this paper, we have studied the delayed secondary emis- 
sion in detail. We have evaluated delayed GRB spectra in 
most detail by numerical simulations, which enable us to 
treat cascade processes including multiple pair creation and 
IC scattering. We have also calculated delayed GRB spectra 
using Eq. (0]i and compared both results. As seen in the pre- 
vious section, both methods agree with each other basically. 

We have especially focused on effects of the CIB. CIB pho- 
tons play a role not only on absorbing high-energy gamma- 
rays but also on being up-scattered as seed photons by cre- 
ated high-energy pairs. USIB photons have larger energy 
than USMB photons, so that USIB photons are more subject 
to duplicated absorption by CIB photons. The USIB com- 
ponent is more sensitive to the CIB than the USMB com- 
ponent. One of the most frequently discussed method to 
probe the CIB is measuring the depletion due to the CIB in 
prompt GRB spectra. However, owing to uncertainties in 
GRB intrinsic spectra, the depletion is hard to be estimated. 
On the other hand, the USIB signature could be more use- 
ful to probe the CIB almost irrespective of intrinsic GRB 
spectra as shown in Fig. 9 for models A and C. While we 
focus on secondary emissions from GRBs, secondary emis- 
sions from blazars are more exp ected (lAharonian et al. 19941 : 
iDai et al. 20021 : iFan et al. 2004h . Our results can be also ap- 
plied to blazars and effects of the CIB can be importrant. 

In this paper, we have clarified the USIB effect especially. 
As shown in Figs. 10-13, not only GLAST but MAGIC and 
VERITAS could detect such high-energy gamma-rays aris- 
ing from USIB photons for sufficiently low redshift bursts. 
There are several characteristic features related to USIB pho- 
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tons: (1) The USIB signature such as high-energy slope ap- 
pears for bursts with the sufficiently low intrinsic high-energy 
cutoff £ , ™ x and/or sufficiently low redshift. (2) If (1) is sat- 
isfied (which means duplicated absorption can be neglected), 
the ratio of the USIB component to the USMB one basically 
reflects the ratio of the CIB energy density to the CMB one. 
(3) The cutoff of delayed spectra is sensitive to the CIB, be- 
cause secondary photons can also be absorbed by CIB pho- 
tons again. Therefore, the energy range, in which USIB pho- 
tons are prominent, becomes narrower, as the source redshift 
becomes higher. (4) As long as the effective intrinsic high- 
energy cutoff £'™ ax is smaller than TeV, the shape of delayed 
secondary spectra, especially the 100 GeV-GeV ratio, is not 
so sensitive to the shape of prompt primary spectra. For ex- 
ample, high-energy USIB slope can be found in both model 
A and model C as seen in the previous section. 

It is important to know the highest energy of gamma- 
rays. These gamma-rays, which may arise from electron syn- 
chrotron radiation, electron SSC, proton synchrotron and par- 
ticles generated by photomeson or photopair production, will 
suffer from absorption due to electron-positron pair creation. 
Numerical simulations in Asano & Inoue (2007) show that the 
intrinsic high-energy cutoff detemined by photon absorption 
is approximated as 

where St is the variability time scale in prompt emission. 
Hence, it is possible to constrain a bulk Lorentz factor by ob- 
serving the highest gamma-ray energy. However, the high- 
est gamma-ray energy can be seen in the prompt spectrum 
only if attenuation due to the CIB is not significant. Because 
delayed secondary emission is also influenced by the high- 
est gamma-ray energy (for example, the typical energy of the 
USMB bump is affected by the intrinsic high-energy cutoff 
Zi™ ax ), observations of these delayed signals could provide us 
with useful information on the source. 

Now, MAGIC has continue d observations and given up- 
per limits for some events (IMAGIC Collaboration 2 006a; 
MAG IC Collaboration 2006bl) . U pper limits are also set by 
other detectors such as Whipple (IConn aughton et al. 1997b 
and STACEE dJarvis et al. 20051) . So far, no excess event 
above ^100 GeV was detected, neither during the prompt 
emission phase nor during the early afterglow. The upper lim- 
its between 85 and 1000 GeV are derived and compatible with 
a naive extension of power-law spectra. 

As we have discussed, the USIB signature in delayed spec- 
tra can be found for GRBs with the relatively low redshift 
and/or relatively low high-energy cutoff. For example, the 
rate of GRBs within z = 0.2-0.3 is ~ afew yr" 1 , so that we 
could see such delayed signals by detectors such as MAGIC 
in the near future. However, the event rate would be not so 
large. The various conditions such as a field of view and ob- 
servational conditions of the detector will reduce the expected 
event rate significantly. Furthermore, the real rate may be 
smaller because the number of bursts that can emit such high- 
energy gamma-rays may be limited. GLAST will see the sig- 
nificant number of GRBs emitting high-energy gamma-rays 
and provide us with information on high-energy spectra of 
prompt emission. In addition, by using GLAST as the moni- 
tor of GRBs that emit high-energy gamma-rays, opportunities 
to observe high energy gamma-rays by ground telescopes will 
also be increased. Since GLAST may see many bursts, some 
of which may include ones with delayed ~ GeV components. 



Furthermore, MAGIC-II is now being constructed and VER- 
ITAS has started observations. These advances in detectors 
could enable us to expect more and more chances to high- 
energy gamma-ray signals from GRBs. 

Here, we discuss possible complications. One is possible 
influence of environments around GRBs. We have neglected 
effects of the environments around host galaxies. But, there 
might be influences from environments such as magnetic field 
of the host galaxies. The second is the existence of the IGM 
field. In this paper, we have assumed the weak IGM field 
with B < 10~ (18-I9) G to estimate delayed secondary fluxes, 
where the magnetic deflection time is not a dominant time 
scale and the angular-size spreading of delayed gamma-rays is 
sufficiently small. Although such a weak magnetic field might 
be possible in the void region, it becomes difficult to observe 
delayed signals when the IGM field is strong enough. This 
is also the reason why the orphan delayed emission (which 
can be expected when 6b > 6j) is difficult to be detected. Es- 
pecially for B > 10~ 16 G, delayed secondary photons will be 
observed as isotropic diffuse signals that are difficult to be 
observed. 

The third possible complication would arise from pro- 
longed intrinsic high-energy emission. TeV signals are also 
expected in the context of the afterglow theory. Although 
the discrimination might be difficult, time-dependent observa- 
tions could enable us to distinguish between two signals, be- 
cause the time evolution will be different between the prompt 
emission and afterglow emission. Other possible late ac- 
tivities discovered by Swift would make further contamina- 
tion. For example, flares can be accompanied not only by 
neutrinos and gamma-rays associated with flares themselves 
dMurase & Nagataki 2006b) but al so by gamma-ray s which 
forward shock electrons up-scatter (Wan g et al. 20061) . 

Finally, we shall comment on high-energy emissions from 
low-lumin osity (LL) GRBs. The recent discovery of XRF 
060218 (ICampana et al. 2 006) implies that there may be a 
different population from usual cosmological high-luminosity 
(HL) GRBs. These LL GRBs are more frequent than usual 
HL GRBs. If true, the high-energy neutrino background 
from LL GRBs can be comparable with tha t from HL GRBs 
dMurase et al. 2006t lGupta & Zhang 20071) . Similarly, we 
can expect the gamma-ray background from LL GRBs is com- 
parable with that from HL GRBs. This is also pointed out by 
Casanova et al. (2007) and Dermer (2006). However, the 
diffuse gamma-ray background from GRBs is much smaller 
than the EGRET bound and the diffuse component will be 
very difficult to be detected. Furthermore, it might be difficult 
to emit TeV photons from LL GRBs, unless they have large 
Lorentz factors. Although we cannot deny the possibility for 
LL GRBs to emit TeV photons so far, we do not consider such 
cases in this paper. 

In summary, we have done the most detailed numerical cal- 
culations, and justified the frequently used approximate ap- 
proach. We have especially studied CIB effects on delayed 
gamma-ray spectra from GRBs in detail. CIB photons are 
important because primary TeV gamma-rays are absorbed by 
them. In addition, we have emphasized that not only CMB 
photons but also CIB photons will be up-scattered by pro- 
duced electron-positron pairs, and enable delayed secondary 
spectra to extend to higher energies. Although this USIB ef- 
fect has not been emphasized in previous studies, it is impor- 
tant since we could obtain additional information on the CIB. 



10 



Murase, Asano, and Nagataki 



We thank T. Kneiske for giving us advice to use their CIB 
spectral model. K.M. thank T. Totani. and K. Ioka for com- 
ments on effects of the CIB. K.M. and S.N. also thank J. Gra- 
not, T. Kamae and M. Teshima for helpful comments. K.M. 
and K.A. thank S. Inoue for helpful discussion. We are grate- 
ful to Z.G. Dai and Y.Z. Fan for useful advices. We would like 
to appreciate for the anonymous referee for giving us prof- 



itable suggestions. 

This work is supported by the Grant-in- Aid for the 21st 
Century COE "Center for Diversity and Universality in 
Physics" from the Ministry of Education, Culture, Sports, Sci- 
ence and Technology (MEXT) of Japan. The work of K.M. is 
supported by a Grant-in- Aid for JSPS Fellows. 



REFERENCES 



Aharonian, F.A. et al. 1994, ApJ, 423, L5 
Aharonian, F.A. et al. 2006, Nature, 440, 1018 
Amenomori, M. et al. 2001, AIPC. Proc., 599, 493 
Ando, S. 2004, MNRAS, 354, 414 
Asano K. 2005, ApJ, 623, 67 
Asano, K., and Inoue, S. 2007, arXiv:0705.2910 
Asano, K., and Nagataki, S. 2006, ApJ, 640, L9 
Asano, K., and Takahara, F. 2003, PASJ, 55, 433 
Atkins, R. 2000, ApJ, 533, LI 19 

Blumenthal, G.R., and Gould, R.J. 1970, Rev. Mod. Phys., 42, 237 

Bruzual, A. G., and Chariot, S. 1993, ApJ, 405, 538 

Campana, et al. 2006, Nature, 442, 1008 

Casanova, S., Dingus, B., and Zhang, B. 2007, ApJ, 656, 306 

Cheng, L.X., and Cheng, K.S. 1996, ApJ, 459, L79 

Connaughton, V. et al. 1997, ApJ, 479, 859 

Dai, Z. G., Zhang, B., Gou, L. J., Meszaros, P., and Waxman, E. 2002, ApJ, 
580, L7 

Dai, Z.G., and Lu, T. 2002, ApJ, 580, 1013 
Dermer, CD. 2006, astro-ph/0610195 

Dermer, C. D., & Atoyan, A. 2003, Phys. Rev. Lett., 91, 071 102 
Dermer, CD., and Atoyan, A. 2004, A&A, 418, L5 
Dermer, CD., Chiang, J., and Mitman, K. 2000, ApJ, 537, 785 
Derishev, E. V., Kocharovsky, V. V., and Kocharovsky, VI. V. 2001, A&A, 
372, 1071 

Enomoto, R. et al. 2002, Astropart. Phys. 16, 235 
Fan, Y. Z., Dai, Z. G, Wei, D. M. 2004, A&A, 415, 483 
Guetta, D., Spada, M., & Waxman, E. 2001, ApJ, 559, 101 
Guetta, D., and Granot, J. 2003, ApJ, 585, 885 

Guetta, D., Hooper, D., Alvarez-Muniz, J., Halzen, E, & Reuveni, E. 2004, 

Astropart. Phys., 20, 429 
Guetta, D., Perna, R., Stella, L., and Vietri, M. 2004, ApJ, 615, L73 
Guetta, D., and Piran, T. 2007, JCAP, 07, 003 
Gupta, N., and Zhang, B. 2007, Astropart. Phys., in press 
Hartmann, D.H. et al. 2002, AIPC. Proc., 662, 477G 
Hauser, M.G, and Dwek, E. 2001, ARA&A, 39, 249 
Hinton, J. A. 2004, New Astron. Rev., 48, 331 
Holder, J. et al. 2006, astro-ph/061 1598 
Hurley, K. et al., 1994, Nature, 372, 652 

Jarvis, B. et al. 2005, Proc. of the 29th International Cosmic Ray Conference, 
4, 455 



Kamae, T. 2006, Private Communication 
Kashlinsky, A. 2006, New Astron. Rev., 50, 208 

Kneiske, T.M., Mannheim, K, and Hartmann, D.H. 2002, A&A, 386, 1 

Kneiske, T.M. et al. 2004, A&A, 413, 807 

Le, T., and Dermer, CD. 2007, ApJ, 661, 394 

Liang, E. et al. 2007, 662, 1111 

Lithwick, Y, and Sari, R. 2001, ApJ, 555, 540 

MAGIC Collaboration, 2006, ApJ, 641 L9 

MAGIC Collaboration, 2006, astro-ph/061 2548 

Meszaros, P. 2006, Rep. Prog. Phys., 69, 2259 

Meszaros, P., Rees, M. J. 1994, MNRAS, 269, L41 

Meszaros, P., Rees, M. J., Papathanassiou, H. 1994, ApJ, 432, 181 

Milagro Collaboration 2007, arXiv:0 705~1554l 

Mirzoyan et al. 2005, Proc. of the 29th International Cosmic Ray Conference, 
4,23 

Murase, K, and Nagataki, S. 2006a, Phys. Rev. D, 73, 063002 
. 2006b, Phys. Rev. Lett., 97, 051101 

Murase, K, Ioka, K., Nagataki, S., and Nakamura, T. 2006, ApJ, 651, L5 

Papathanassiou, H., and Meszaros, P., ApJ, 471, L91 

Peer, A., and Waxman, E. 2004, ApJ, 613, 448 

Perry, D. et al. 1999, A&A, 138, 601 

Poirier, J. at al. 2003, Phys. Rev. D, 67 

Porciani, C, and Madau, P. 2001, ApJ, 548, 522 

Plaga, R. 1995, Nature, 374, 30 

Razzaque, S., Meszaros, P., and Zhang, B. 2004, ApJ, 613, 1072 
Sreekumar, P. et al. 1998, ApJ, 494, 523 

Stecker, F.W., de Jager, O.C., and Salamon, M.H. 1992, ApJ, 390, L49 
Stecker, F.W., Malkan, M.A., and Scully, S.T. 2006, ApJ, 648, 774 
Sommer, M. 1994, ApJ, 422, L63 

Strong, A.W., Moskalenko, I.V., and Reimer, 0. 2004, ApJ, 648, 774 

Totani, T. 1998, ApJ, 509, L81 

Totani, T., and Takeuchi, T.T. 2002, ApJ, 570, 470 

Wang, X.Y., Dai, Z.G., and Lu, T. 2001, ApJ, 556, 1010 

Wang, X.Y., Cheng, K.S., Dai, Z.G., and Lu, T. 2004, ApJ, 604, 306 

Wang, X.Y., Li, Z., and Meszaros, P. 2006, ApJ, 641, L89 

Waxman, E. and Bahcall, J. 1997, Phys. Rev. Lett., 78, 2292 

Zhang, B. 2007, ChJAA, 7, 1 

Zhang, B., and Meszaros, P. 2001, ApJ, 559, 110 



